EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 




CERN-PH-EP-20 12-368 
18 Dec 2012 



Centrality determination of Pb-Pb collisions 
at VsSn= 2 - 7 6 TeV with ALICE 

CO ALICE Collaboration 

i— i 

O 

<N 

Ch Abstract 

qq This publication describes the methods used to measure the centrality of inelastic Pb-Pb 

i— l collisions at a center-of-mass energy of 2.76 TeV per colliding nucleon pair with ALICE. 

The centrality is a key parameter in the study of the properties of QCD matter at extreme 
temperature and energy density, because it is directly related to the initial overlap region 
V of the colliding nuclei. Geometrical properties of the collision, such as the number of 

q participating nucleons and number of binary nucleon-nucleon collisions, are deduced from 

a Glauber model with a sharp impact parameter selection, and shown to be consistent with 
i those extracted from the data. The centrality determination provides a tool to compare 

ALICE measurements with those of other experiments and with theoretical calculations. 
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1 Introduction 

Ultra-relativistic heavy-ion collisions at the Large Hadron Collider (LHC) produce strongly 
interacting matter under extreme conditions of temperature and energy density, similar to those 
prevailing in the first few microseconds after the Big Bang (U. 

Since nuclei are extended objects, the volume of the interacting region depends on the impact 
parameter (b) of the collision, defined as the distance between the centers of the two colliding 
nuclei in a plane transverse to the beam axis. It is customary in the field of heavy-ion physics 
to introduce the concept of the centrality of the collision, which is directly related to the impact 
parameter, and infered by comparison of data with simulations of the collisions. 

The purely geometrical Glauber model [2], which typically is used in this context, has its ori- 
gins in the quantum mechanical model for p-A and A-A scattering described in 0]-[5]]. The 
model treats a nuclear collision as a superposition of binary nucleon-nucleon interactions. The 
volume of the initial overlap region is expressed via the number of participant nucleons. A 
participant nucleon of one nucleus is defined as a nucleon that undergoes one or more binary 
collisions with nucleons of the other nucleus. The number of participants and spectators, 7V par t 
and A^ S pec = 2A — 7V part , where A is the total number of nucleons in the nucleus (mass number), 
and the number of binary collisions N co \\ are calculated for a given value of the impact parameter 
and for a realistic initial distribution of nucleons inside the nucleus, and assuming that nucleons 
follow straight trajectories. This approach provides a consistent description of p-A, d-A, and 
A-A collisions, and is especially useful when comparing data from different experiments or 
from different collision systems and to theoretical calculations. 

Neither the impact parameter nor geometrical quantities, such as Af part , Af spec , or N co \\ are directly 
measurable. Two experimental observables related to the collision geometry are the average 
charged-particle multiplicity 7V c h and the energy carried by particles close to the beam direc- 
tion and deposited in Zero-Degree Calorimeters (ZDC), called the zero-degree energy Ezdc- 
The average charged-particle multiplicity is assumed to decrease monotonically with increas- 
ing impact parameter. The energy deposited in the zero-degree calorimeters, Ezdc, is directly 
related to the number of spectator nucleons Af spec , which constitute the part of the nuclear vol- 
ume not involved in the interaction. However, unlike 7V c h, Ezdc does not depend monotonically 
on the impact parameter b because nucleons bound in nuclear fragments with similar magnetic 
rigidity as the beam nuclei remain inside the beam-pipe and therefore are not detected in the 
ZDC. Since fragment formation is more important in peripheral collisions, the monotonic re- 
lationship between Ezdc and b is valid only for relatively central events (small b). For this 
reason, the zero-degree energy measurement needs to be combined with another observable 
that is monotonically correlated with b. 

The centrality is usually expressed as a percentage of the total nuclear interaction cross section 
a 0. The centrality percentile c of an A-A collision with an impact parameter b is defined by 
integrating the impact parameter distribution do/db as 



In ALICE, the centrality is defined as the percentile of the hadronic cross section corresponding 
to a particle multiplicity above a given threshold (N^ R ) or an energy deposited in the ZDC 
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below a given value (E^c) in the ZDC energy distribution do/dE : 
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The procedure can be simplified by replacing the cross section with the number of observed 
events, corrected for the trigger efficiency. However, at LHC energies, the strong electromag- 
netic fields generated by the heavy ions moving at relativistic velocity lead to large cross sec- 
tions for QED processes (6H91. Although the cross sections for these processes exceed those for 
the hadronic cross section by several oders of magnitude, they only contaminate the hadronic 
cross section in the most peripheral collisions. For this reason one may choose to restrict the 
centrality determination to the region where such contamination is negligible. The fraction of 
hadronic events excluded by such cut as well as the trigger efficiency can be estimated using a 
model of the nuclear collision and the related particle production. 

In this paper, we report on the centrality determination used in the analyses of the Pb-Pb col- 
lision data from the 2010 and 2011 run recorded with the ALICE detector ifTOll . Specifically 
the analysis presented here is done with a subset of the 2010 data, but the methods and results 
are valid for 2011 as well. In Section |2j we describe the implementation of the Glauber model 
used by ALICE. We extract mean numbers of the relevant geometrical quantities for typical 
centrality classes defined by classifying the events according to their impact parameter. Sec- 
tion [3] describes the experimental conditions and the event selection with particular emphasis 
on the rejection of QED and machine-induced backgrounds. Section [4] presents the methods 
employed by ALICE for the determination of the hadronic cross section, needed for the abso- 
lute determination of the centrality. The main method uses the VZERO amplitude distribution 
fitted with the Glauber model. The systematic uncertainty is obtained by comparing the fit to 
an unbiased VZERO distribution obtained by correcting the measured one by the efficiency of 
the event selection and the purity of the event sample. Section [5] presents the determination 
of the centrality classes using either the multiplicity at mid-rapidity or the energy deposited in 
the ZDC. We discuss the relation between the measured multiplicity and geometrical quanti- 
ties connected to centrality, established by the Glauber Model. These are nearly identical to 
those obtained in Section [2} classifying the events according to their impact parameter, which 
are therefore used as reference in all ALICE analyses. Section [6] presents the precision of the 
centrality determination in ALICE. Section [7] summarizes and concludes the paper. 



The Glauber model is widely used to describe the dependence of A^ art and N co \\ on b in p-A, d-A 
and A-A collisions [EHSJ. The purpose of Monte Carlo implementations of the Glauber model 
lfT9l l20l is to compose two nuclei out of nucleons and simulate their collision process event- 
by-event. Geometrical quantities are calculated by simulating many nucleus-nucleus collisions. 
Mean values of these quantities are calculated for centrality classes defined by classifying the 
events according to their impact parameter b. 

Following ED, the first step in the Glauber Monte Carlo is to prepare a model of the two nuclei 
by defining stochastically the position of the nucleons in each nucleus. The nucleon position 
in the 208 Pb nucleus is determined by the nuclear density function, modeled by the functional 
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Fig. 1: Compilation of total a^, elastic <7^ N , and inelastic c^n cross sections of pp and pp collisions 

lfTTlfl3ll. The o£' N curve is a fit performed by the COMPETE Collaboration also available at EICH). The 

pp data from ATLAS [151, CMS ESI, TOTEM, 021 and ALICE HI agree well with the interpolation 
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form (modified Woods-Saxon or 2-parameter Fermi distribution): 



p(r) 



l+w(r/R) 2 
P °l+exp(^) 



(3) 



The parameters are based on data from low energy electron-nucleus scattering experiments ll22l . 
Protons and neutrons are assumed to have the same nuclear profile. The parameter po is the 
nucleon density, which provides the overall normalization, not relevant for the Monte Carlo 
simulation, R = (6.62 ± 0.06) fm is the radius parameter of the 208 Pb nucleus and a = (0.546 ± 
0.010) fm is the skin thickness of the nucleus, which indicates how quickly the nuclear density 
falls off near the edge of the nucleus. The additional parameter w is needed to describe nuclei 
whose maximum density is reached at radii r > (w = for Pb). In the Monte Carlo procedure 
the radial coordinate of a nucleon is randomly drawn from the distribution 4nr 2 p(r) and po 
is determined by the overall normalization condition / p(r)d 3 r = A. We require a hard-sphere 
exclusion distance of d m [ n = 0.4 fm between the centers of the nucleons, i.e. no pair of nucleons 
inside the nucleus has a distance less than d m [ n . The hard- sphere exclusion distance, character- 
istic of the length of the repulsive nucleon-nucleon force, is not known experimentally and thus 
is varied by 100% (d min = (0.4 ±0.4) fm). 

The second step is to simulate a nuclear collision. The impact parameter b is randomly selected 
from the geometrical distribution dP/db ~ b up to a maximum Z? max ~ 20 fm > 27?pt,. The 
maximum value of impact parameter £> max is chosen large enough to simulate collisions until 
the interaction probability becomes zero. This is particularly important for the calculation of the 
total Pb-Pb cross section. The nucleus-nucleus collision is treated as a sequence of independent 
binary nucleon-nucleon collisions, where the nucleons travel on straight-line trajectories and the 
inelastic nucleon-nucleon cross section is assumed to be independent of the number of collisions 
a nucleon underwent previously, i.e. the same cross section is used for all successive collisions. 
Two nucleons from different nuclei are assumed to collide if the relative transverse distance 
between centers is less than the distance corresponding to the inelastic nucleon-nucleon cross 
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section d < y Ojj^j 1 / 7T. A Gaussian overlap function can be used as an alternative to the black- 
disk nucleon-nucleon overlap function (23). It makes no significant difference within systematic 
uncertainty in the global event properties. 



Table 1: Geometric properties (Af par t ? Af co n, 7aa) of Pb-Pb collisions for centrality classes defined by 
sharp cuts in the impact parameter b (in fm). The mean values, the RMS, and the systematic uncertainties 
are obtained with a Glauber Monte Carlo calculation. 
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The number of collisions 7V co ii and the number of participants Af part are determined by count- 
ing, respectively, the binary nucleon collisions and the nucleons that experience at least one 
collision. Following the notation in 0, the geometric nuclear overlap function Taa is then 
calculated as Taa — N co \\/ cr^N' anc ^ represents the effective nucleon luminosity in the collision 
process. 

For nuclear collisions at y^NN = 2.76 TeV, we use O"^ 1 = (64 ±5) mb, estimated by inter- 
polation [11] of pp data at different center-of-mass energies and from cosmic rays lfT2l [T4l . 
and subtracting the elastic scattering cross section from the total cross section. The interpo- 
lation is in good agreement with the ALICE measurement of the pp inelastic cross section at 
= 2.76 TeV, cr^ 1 = (62.8 ±2.4+J J) mb d, and with the measurements of ATLAS 
Gl, CMS IH, and TOTEM G3 at ^nn= 7 TeV, as shown in Fig.|T] 
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The total Pb-Pb cross section is calculated as OpbPb = N eyt (N C0 \\ > l)/N &v t(N co vi. > 0) x 7rZ?^ ax , 
i.e. the geometrical value corrected by the fraction of events with at least one nucleon-nucleon 
collision. We obtain OpbPb = (7 .64±0.22(syst .)) b, in agreement with the ALICE measurement 
dPbPb = (lJ±OA(stat.)t° j(syst.)) b & 

Table [T] report the mean number of participants (Af part )and collisions (N co \\), and the mean nu- 
clear thickness function (7aa) f° r centrality classes defined by sharp cuts in the impact pa- 
rameter b calculated with the Glauber model (Fig. [2]). The root mean square (RMS) of these 
distributions is a measure for the magnitude of the dispersion of the quantities. 

The systematic uncertainties on the mean values are obtained by independently varying the 
parameters of the Glauber model within their estimated uncertainties. More specifically, the 
default value of the nucleon-nucleon cross section of cr^ 1 = 64 mb was varied between 59 mb 
and 69 mb. The Woods-Saxon parameters were varied by one standard deviation to determine 
uncertainties related to the nuclear density profile. The minimum distance of 0.4 fm between 
two nucleons of the same nucleus was varied by 100%, from to 0.8 fm to evaluate the effects 
of a nucleon hard core (as mentioned above). Figure [3] shows the resulting variations for Pb-Pb 
collisions at v /Snn= 2.76 TeV. The total systematic uncertainty reported in Table[T]was obtained 
by adding in quadrature the deviations from the default result for each of the variations listed 
above. The uncertainty of A^ art ranges from about 3-4% in peripheral collisions to < 1% in 
central collisions, the uncertainty of N co \\ ranges from about 7% in peripheral collisions to about 
11% in central collisions, the uncertainty of Taa ranges from about 6% in peripheral collisions 
to about 3% in central collisions. The nuclear overlap function Taa is often used to compare 
observables related to hard processes in A-A and pp collisions. Since Taa — N co \\/o^, it has 
the same systematic uncertainties as N co \\ except that the uncertainty on o^ 1 cancels out. 

Finally it is worth to note that more sophisticated implementations of the Glauber model ll23l - 
l25ll suggest that effects not included in our Glauber model, such as the changes of the excluded 
volume on the nuclear density and two-body correlations, can be approximated by slightly 
adjusting the Woods-Saxon parameters. The modified parameters, however, are well covered 
by the systematic uncertainty quoted above for the parameters that we use. 




b (fm) N Part 



Fig. 2: Geometric properties of Pb-Pb collisions at y^NN = 2.76 TeV obtained from a Glauber Monte 
Carlo calculation: Impact parameter distribution (left), sliced for percentiles of the hadronic cross sec- 
tion, and distributions of the number of participants (right) for the corresponding centrality classes. 
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Fig. 3: Sensitivity of A^ part (left) and A^oii (right) to variations of parameters in the Glauber Monte Carlo 
model of Pb-Pb collisions at - v /5 , nn= 2.76 TeV. The gray band represents the RMS of Af part and N co \\ 
respectively It is scaled by a factor 0.1 for visibility. 

3 Experimental Conditions 
3.1 The ALICE detector 

ALICE is an experiment dedicated to the study of heavy-ion collisions at the LHC. A detailed 
description of the apparatus is given in Ref. [10J. Here, we briefly describe the detector compo- 
nents used in this analysis. 

The Silicon Pixel Detector (SPD) is the innermost part of the Inner Tracking System (ITS). 
It consists of two cylindrical layers of hybrid silicon pixel assemblies positioned at an average 
radial distances of 3.9 and 7.6 cm from the beam line, with a total of 9.8 x 10 6 pixels of size 50 x 
425 jUm 2 , read out by 1200 electronic chips. The SPD coverage for particles originating from 
the center of the detector is \rf | < 2.0 and |rj | < 1.4 for the inner and outer layers, respectively. 
Each chip provides a fast signal if at least one of its pixels is hit. The signals from the 1200 
chips are combined in a programmable logic unit which supplies a trigger signal. The fraction 
of SPD channels active during 2010 data taking was 70% for the inner and 78% for the outer 
layers. 

The VZERO detector consists of two arrays of 32 scintillator cells placed at distances z = 
3.4 m and z = —0.9 m from the nominal interaction point, along the beam line, covering the 
full azimuth. The VZERO detector is within 2.8 < r/ < 5.1 (VZERO-A) and -3.7 < i] < 
— 1.7 (VZERO-C). Both amplitude and time of signals in each scintillator are recorded. The 
VZERO time resolution is better than 1 ns, allowing discrimination of beam-beam collisions 
from background events produced upstream of the experiment. The VZERO is also used to 
provide a trigger signal (see |3.2| ). 

The Time Projection Chamber (TPC) is used for charged particle trajectory reconstruction, track 
momentum measurement and particle identification. The ALICE TPC is a large cylindrical drift 
detector whose active volume extends radially from 85 to 247 cm, and from -250 to +250 cm 
along the beam direction. The active volume of nearly 90 m 3 is filled with a gas mixture of Ne 
(85.7%), C0 2 (9.5%) and N 2 (4.8%) until the end of 2010, and Ne (90%) and C0 2 (10%) since 
the beginning of 2011. A central electrode maintained at -100 kV divides the TPC into two 
sections. The end-caps are equipped with multiwire proportional chambers with cathode pad 
readout. For a particle traversing the TPC up to 159 position signals (clusters) are recorded. 
The cluster data are used to reconstruct the charged particle trajectory as well as to calculate the 
particle's specific energy loss used to identify the species of the particle which has produced the 
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track. 

The two Zero Degree Calorimeters (ZDC) in the ALICE experiment measure the energy of 
spectator (non-interacting) nucleons: ZP measures protons and ZN measures neutrons. They 
are situated about 114 m from the interaction point on each side of the experiment [fTOl . Each 
ZDC consists of two quartz fiber sampling calorimeters: the neutron calorimeter positioned 
between the two beam pipes downstream of the first machine dipole that separates the two 
charged particle beams, and the proton calorimeter positioned externally to the outgoing beam 
pipe. The energy resolution at beam energy is estimated to be 20% for the neutron (20.0% for 
ZNC, 21.2% for ZNA) and 24% for the proton calorimeters, respectively. 
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Fig. 4: Time distribution of signals in the VZERO detector on the A side. The peaks corresponding to 
beam-beam, beam-gas and satellite collision events are clearly visible. 

3.2 Data set and online event selection 

During the first LHC Pb-Pb run in 2010, beams of four bunches with about 10 7 Pb ions per 
bunch collided at v /snn = 2.76 TeV, with an estimated luminosity of 5 x 10 23 cm _2 s _1 . AL- 
ICE collected about 90 million nuclear collision events using different interaction triggers with 
increasingly tighter conditions. These triggers used VZERO and SPD detector signals in coin- 
cidence with a bunch crossing corresponding to a beam-beam collision: 

— V0AND: signals in VZERO-A and VZERO-C; 

— 3-out-of-3: signals in VZERO-A and VZERO-C and at least 2 chips hit in the outer layer 
of the SPD; 

— 2-out-of-3\ two of the three conditions listed above. 

The threshold in the VZERO detector for each of the VZERO tiles corresponded approximately 
to the energy deposition of one minimum ionizing particle. 

Control events were also collected with the same trigger logic, in coincidence with only one 
beam crossing the ALICE interaction point (from either the A or the C side) or with no beam at 
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all ("empty"). The luminous region had an RMS width of 5.9 cm in the longitudinal direction 
and 50 /im in the transverse direction. For the estimated luminosity, using the least selective of 
the interaction triggers, the observed rate was about 50 Hz. This was mainly due to electromag- 
netically induced processes [26] . These processes have large cross sections at LHC energies 
but generate low multiplicities and therefore do not contribute to the typical particle multiplici- 
ties of interest for the present paper. The trigger rate without beam was negligible and the rate 
in coincidence with bunches of only one beam was about 1 Hz. The probability for collision 
pile-up per triggered event was less than 10 -4 . 




Fig. 5: Correlation between the sum and the difference of times recorded by the neutron ZDC on either 
side of the interaction region. The large cluster in the middle corresponds to collisions between ions in 
the nominal RF buckets of each beam, while the small clusters along the diagonals (spaced by 2.5 ns in 
the time difference) correspond to collisions in which one of the ions is displaced by one or more RF 
buckets. 



3.3 Offline event selection 

The offline event selection is applied with the purpose of selecting hadronic interactions with the 
highest possible efficiency, while rejecting the machine-induced and physical backgrounds. The 
offline event selection replays the on-line trigger condition, using the same quantities calculated 
offline, so that events triggered by noise in the SPD are discarded, and the weighted time average 
over all channels is used for the VZERO, leading to a better time resolution. In addition, the 
offline event selection rejects the machine-induced background and parasitic collisions. This 
contamination amounts to about 25% of all collected events. To keep the conditions of all 
detectors as uniform as possible (in particular those around mid-rapidity, such as the SPD), the 
centrality analysis was restricted to a region around the vertex, \z v tx\ % 10 cm. 

3.3.1 Machine-Induced Background 

One source of machine-induced background is due to beam-gas events, caused by one of the 
beam interacting with the residual gas in the beam-pipe; another source of background are 
events where ions in the beam halo interact with mechanical structures in the machine. These 
interactions mostly occur outside of the interaction region and thus produce a signal that is "too 
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early" in the same-side VZERO, compared to a collision that occurs in the nominal interaction 
region between the VZERO detectors. Therefore these events can be rejected using the timing 
information of the VZERO. This is illustrated in Fig. [4] which shows the arrival time of particles 
at the VZERO A detector relative to the nominal beam crossing time. Beam-halo or beam-gas 
interactions are visible as secondary peaks in the time distribution because particles produced 
in background interactions arrive at earlier times in the detector relative to particles produced in 
beam-beam collisions at the nominal vertex, which are the majority of the signals. Other small 
peaks between these main ones arise from satellite collisions. 

Another source of machine-induced background is due to parasitic collisions from debunched 
ions. The radio-frequency (RF) structure of the LHC of 400 MHz is such that there are 10 
equidistant RF buckets within the 25 ns time interval between two possible nominal bunch 
positions. Therefore the buckets are spaced by 2.5 ns. Only one of them should be populated by 
ions [27 j. However, ions can "jump" into one of the neighboring buckets. Therefore collisions 
occur either between ions in the nominal RF buckets but also between one or two ions displaced 
by one or more RF buckets. This causes a displacement in the Z- vertex position of 2.5 ns/2 • c — 
37.5 cm, well outside the fiducial region \z v tx\ ^5 10 cm. Those events are thus to be considered 
as "background" and are rejected using the correlation between the sum and the difference of 
times measured in each of the neutron ZDCs, as shown in Fig. [5j Such satellite collisions can 
also be rejected using the vertex cut. 

After the event selection, the remaining machine-induced background, estimated from the con- 
trol triggers (i.e. triggers that fire for coincidences between empty and filled or empty and empty 
bunches), is negligible. 




Ezna (TeV) 

Fig. 6: Correlation between signals in the two neutron zero-degree calorimeters, ZNA and ZNC. The 
figure is taken from (9J. Single electromagnetic dissociation events produce signal in only one of the 
calorimeters. Mutual dissociation and hadronic interactions populate interior of the plot and can be 
distinguished from each other by the signal in ZEM. 
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3.3.2 Electromagnetic interaction background 

At the LHC energy, the cross sections for electromagnetic (EM) processes, generated by the 
EM fields of relativistic heavy ions, are enormous (^(kbarn)) (6l-[9l. This is the main physical 
background, and needs to be rejected in heavy-ion collisions to isolate hadronic interactions. 
QED processes consist of: photo-production and photo-nuclear interactions. Photo-production 
results in the creation of an e + e~ pair. Photo-nuclear interactions, where one photon from the 
EM field of one of the nuclei interacts with the other nucleus, possibly fluctuating to a vector 
meson, yield a low multiplicity of soft particles in the ALICE central barrel. In the case of 
single photo-production the particle multiplicity is asymmetric within the event. Along the 
beam direction the electromagnetic dissociation (EMD) cross sections are large resulting in a 
non-negligible probability for one neutron emission from either nucleus. 

The EMD cross sections have been measured in a special run triggering on a signal in one of 
the neutron ZDCs, ZNA or ZNC, with a threshold placed well below the single neutron signal 
to detect the neutrons from Giant Dipole Resonance (GDR) decay emitted very close to beam 
rapidity (21. The recorded event sample is dominated by electromagnetic dissociation of one 
or both nuclei measured to be <j sin g leEMD = 187.4 ±02(stat.)t\l'l(syst.) b compared to the 
mutual EMD cross section of a mutualEMD = 5 J ±0.1 (st at.) ±0A(syst.) b. The single EMD 
events can be clearly identified when correlating the response of ZNA and ZNCs (Fig. [6]). The 
additional requirement of a signal in an electromagnetic calorimeter close to beam rapidity 
(ZEM) allows one to distinguish between mutual EMD and hadronic interaction events. 

In order to reduce the contribution due to single neutron emission, we require a ZDC signal 
three standard deviations above the single neutron peak. This selection rejects about 3% of 
all events within 10 cm from the nominal interaction point after removal of beam-gas and par- 
asitic collisions, and only removes events for peripheral collisions (in the 90-100% region). 
The coincidence of the ZDC signals rejects all the single neutron emission events. The simul- 
tanous emissions from both nuclei still are accepted, which, however, are only relevant for very 
peripheral collisions. 

For systematic studies, another selection based on the information from the TPC is used, where 
at least one track reconstructed in the TPC is requested in order to keep the event. This selection 
removes few peripheral hadronic interactions, and strongly suppresses the EM background. 



4 Determination of the hadronic cross section 

In order to classify the collisions in percentiles of the hadronic cross section using the charged 
particle multiplicity, it is necessary to know the particle multiplicity at which the purity of the 
event sample and the efficiency of the event selection becomes 100%. We define the Anchor 
Point (AP) as the amplitude of the VZERO detector equivalent to 90% of the hadronic cross 
section, which determines the absolute scale of the centrality. The determination of the AP 
requires the knowledge of the trigger efficiency and the remaining background contamination 
in nuclear collision events. Two methods have been used to study this. The difference in the 
results obtained with the two methods is used to estimate the systematic uncertainty by 



Simulating the multiplicity distribution (see Sec. |4.1| ). In the first approach, we use a full sim- 
ulation of hadronic and EM processes, including a detailed description of the detector response, 
to study the efficiency of the event selection (Sec. |4.1.1[ ) and to estimate the background con- 
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Fig. 7: Efficiency of the three online triggers (2— out — of— 3, V0AND, 3— out — of— 3) used for Pb-Pb 
collisions as a function of the VZERO amplitude calculated with HIJING and AMPT, and measured in 
dedicated pp runs. The efficiency in the simulation has been calculated for events with A^ part = 2. 



tamination (Sec. |4.1.2| ). The real multiplicity distribution, corrected for efficiency and purity, 
allows direct access to the AP. 

Fitting the multiplicity distribution (see Sec. |4.2| ). In the second method, we use the Glauber 
Monte Carlo, combined with a simple model for particle production, to simulate a multiplic- 
ity distribution which is then compared to the experimental one. The simulated distribution 
describes the experimental one down to the most peripheral events where they start to deviate 
due to background contamination and limited trigger efficiency. The location of the divergence 
between the data and simulation can be used to define the AP. 

The centrality determination is performed for different trigger and detector settings. The differ- 
ent triggers change the fraction of accepted events from EM processes, which determines the 
shape of the multiplicity distribution for very peripheral collisions below the AP. The position 
of the AP is very stable for the entire 2010 (and 201 1) run period and does not change within the 
quoted systematic uncertainty discussed below. Small variations in detector conditions induce 
small changes in the position of the edge of the multiplicity distribution for most central events. 
Nevertheless, the centrality determination, adjusted to account for small changes in the detector 
configuration, provides a stable centrality selection for the entire data-taking period (the mean 
fraction of events in the 0-1% bin is 0.01 with a RMS of 0.001, and in the 40-50% bin the mean 
fraction is 0.101 with an RMS of 0.002). 

4.1 Method 1: Correcting the multiplicity distribution 

With this method the AP is determined by evaluating the efficiency of the event selection and 
by estimating the purity of the obtained event sample. 
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4.1.1 Efficiency of the event selection 

The efficiency for the different event selections is studied with simulations of hadronic reactions 
and with dedicated pp runs. For simulations we use HIJING ll28ll or AMPT ll29l with a full 
GEANT [30J description of the ALICE detector and a trigger emulator. In the simulations the 
efficiency is defined as the ratio of events selected by a given condition to all generated events. 
In the two dedicated pp runs, which were taken at the end of the 2010 run at ^fs = 7 TeV, the 
detector conditions were similar as those in the Pb-Pb run. The VZERO gain was adjusted 
such that the response to minimum ionizing particles (MIP) corresponded to 3 ADC channels 
for one run, or 6 for the other run. For the Pb-Pb run, it was set to 4 channels, i.e. between 
the two tested conditions. In the special pp runs, we used a minimum interaction trigger, which 
requires a logical OR between a hit in the SPD and in either of the two VZERO detectors 
(CINT1 trigger condition). The same event selection criteria as used in the Pb-Pb run as the 
trigger have been applied. The relative event selection efficiency is defined as the ratio of events 
selected by a given condition to all the events recorded with the pp minimum-bias interaction 
trigger (CINT1). Since the pp minimum-bias interaction trigger (CINT1) has an efficiency that 
is effectively 100% for non-diffractive events [18], the relative efficiency measured in the pp 
runs, shown in Fig. [7} can be qualitatively compared to that obtained in Pb-Pb simulations with 
Afpart = 2. Except for very low amplitudes, results from HIJING and AMPT are in very good 
agreement. AMPT predicts a slightly higher efficiency (about 0.5%), as a consequence of the 
broader rapidity distribution. The comparison with the pp runs shows a reasonable agreement 
for the "MIP = 6" case, while the "MIP = 3" is clearly lower. 

For the Pb-Pb run, the efficiency of the event selection is calculated using the average of results 
obtained with HIJING and AMPT. The efficiency of the interaction triggers is 99.4%, 97.1%, 
96.9% respectively for 2-out-of-3, V0AND, 3-out-of-3 using HIJING and 99.7%, 98.6%, 98.4% 
using AMPT. The line in Fig. [7} corresponding to the 90% of the hadronic cross section, shows 
that the trigger is always fully efficient for the 90% most central collisions, except for the "MIP 
= 3" pp case, where the efficiency is 95%. 

4.1.2 Remaining contamination 

The purity of the data sample passing a given event selection is estimated using HIJING sim- 
ulations Il28ll for hadronic processes and QED [8] and STARLIGHT Q for the simulations of 
the EM background. For the electromagnetic dissociation we assume that the selection based 
on the signal 3<7 above the single-neutron peak in the ZDCs (see |3.3.2| ) is fully efficient. 

In Fig. [8} data taken with the V0AND interaction trigger are compared to the sum of HIJING and 
background (QED + STARLIGHT) simulations with the same event selection. The simulations 
are scaled to the known cross sections: 

- HIJING (hadronic): a H = 7.66 b ESI; 

- QED (EM): o Q = 92 kb GQ; 

- STARLIGHT (single neutron dissociation): o S ns = 24.2 b Q; 

- STARLIGHT (double neutron dissociation): a SND = 240 mb Q. 

The sum of the simulations is normalized to the data in the region 150 < VZERO amplitude < 
500, where there is no background contamination. The contribution from QED is completely 
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Fig. 8: VZERO amplitude distribution in data (red points) and simulations with the VOAND interaction 
trigger. The data are compared to the sum of HIJING + QED + STARLIGHT simulations (histogram) 
with the same event selection. 

removed by the VOAND trigger. The dashed lines, indicating respectively 80% and 90% of the 
hadronic cross section, show that there is no significant background contamination for colli- 
sions more central than 90%. The region 90-100% is reasonably understood as the agreement 
between data and simulation is quite good. The remaining discrepancy between the data and 
the sum of all contributions is included in the systematic uncertainty. 

To assign a systematic uncertainty, the comparison is made for the three online interaction 
triggers and other event selections requiring (i) VOAND + TPC: one track fully reconstructed in 
the TPC on top of the VOAND trigger; (ii) VOAND + ZDC: 3<7 cut above single neutron peak 
in ZDC on top of the VOAND trigger. For all these variations a cross section is calculated and 
the difference is included in the systematic uncertainty. 

Figure [9] shows the purity of the various Pb-Pb event samples after those selections. The pu- 
rity, plotted as a function of the VZERO amplitude (V), is defined as the fraction of hadronic 
collisions over all the events selected with a given condition: 

dN i Qh 

purity = dYlHNH (4) 

dV \ H N H ^ dV \ sns N S ns dV \ SND N SND ^ dV \Qn q 

where o x and N x are the cross sections and number of events for a given process, x 9 where 
x = H, SNS, SND, and g, for HIJING, STARLIGHT single, STARLIGHT double, and QED, 
respectively. 

The purity of the event sample can be verified using the correlation of the energy deposition in 
the two sides of the ZN calorimeter, similar to the one shown in Figj6} Single neutron peaks are 
visible in the 80-90% centrality class, which may indicate some remaining contamination from 
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Fig. 9: Purity of the 3 online interaction triggers (2-out-of-3, VOAND, 3-out-of-3) and other event se- 
lections used for Pb-Pb collisions as a function of the VZERO amplitude calculated with HIJING, 
STARLIGHT and QED simulations. The dashed line indicates 90% of the hadronic cross section. 

EMD events. However their origin can be also attributed to asymmetric Pb-Pb events, as well 
as a pile up of an EMD and a hadronic collision. Since this contamination can not be easily 
removed, analyses that use peripheral classes like 80-90% assign an additional 6% systematic 
uncertainty on the event selection to take into account the possible contamination from EMD. 
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Fig. 10: Distribution of the sum of amplitudes in the VZERO scintillators. The distribution is fitted with 
the NBD-Glauber fit (explained in the text) shown as a line. The centrality classes used in the analysis 
are indicated in the figure. The inset shows a zoom of the most peripheral region. 



4.2 Method 2: Fitting the multiplicity distribution 

Another independent way to define the AP uses a phenomenological approach based on the 
Glauber Monte Carlo to fit the experimental multiplicity distribution. The Glauber Monte Carlo 
uses the assumptions mentioned above plus a convolution of a model for particle production, 
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based on a negative binomial distribution (NBD). This latter assumption is motivated by the 
fact that in minimum bias pp and pp collisions at high energy, the charged particle multi- 
plicity do/dN c h has been measured over a wide range of rapidity and is well described by a 
NBD OH [32]l. This approach allows one to simulate an experimental multiplicity distribution 
(e.g. VZERO amplitude), which can be compared with the one from data. 



Figure 10 shows the distribution of VZERO amplitudes for all events triggered with the 3-out-of- 



3 trigger (see |3.2| ) after removing the beam background (see |3.3.1| ), part of the EM background 



with the ZDC cut (see 3.3.2 ) and a Z-vertex cut \z v tx\ < 10 cm. The multiplicity distribution has 



the classical shape of a peak corresponding to most peripheral collisions (contaminated by EM 
background and by missing events due to the trigger inefficiency), a plateau of the intermediate 
region and an edge for the central collisions, which is sensitive to the intrinsic fluctuations of 
Npart and d/V c h/drj and to detector acceptance and resolution. 

The Glauber Monte Carlo defines, for an event with a given impact parameter b, the correspond- 
ing Afp ar t and A/coii. The particle multiplicity per nucleon-nucleon collision is parametrized by 
a NBD. To apply this model to any collision with a given 7V part and N co \\ value we introduce 
the concept of "ancestors", i.e. independently emitting sources of particles. We assume that 
the number of ancestors Af ances tors can be parameterized by Af ances tors = / ■ Af par t + (1 — /) • N co \\. 
This is inspired by two-component models (33l [34]], which decompose nucleus-nucleus col- 
lisions into soft and hard interactions, where the soft interactions produce particles with an 
average multiplicity proportional to Af part , and the probability for hard interactions to occur is 
proportional to A^oll- We discuss the independence of the fit results of this assumption below 
(Sec. gZT) . 

To generate the number of particles produced per interaction, we use the negative binomial 
distribution 

P M= nn + k) (ii/k)» 

r(n+i)r(*) " (mA+i)^ 5 

which gives the probability of measuring n hits per ancestor, where ji is the mean multiplicity 
per ancestor and k controls the width. For every Glauber Monte Carlo event, the NBD is sam- 
pled Ancestors times to obtain the averaged simulated VZERO amplitude for this event, which 
is proportional to the number of particles hitting the hodoscopes. The VZERO amplitude dis- 
tribution is simulated for an ensemble of events and for various values of the NBD parameters 
fi, k, and the Ancestors parameter /. A minimization procedure is applied to find the parame- 



ters which result in the smallest # 2 , also shown in Fig. [101 The fit is performed for VZERO 
amplitudes large enough so that the purity of the event sample and the efficiency of the event 
selection is 100%. That leaves a very broad range in the amplitude values that can be fitted to 
extract parameters /, ji and k directly from the data. The amplitude, above which we have 90% 
of the hadronic cross section, defines the AR The quality of the fit is good, as the # 2 /NDF is 
approximately unity for all fits. We note that the high multiplicity tail, which is quite sensitive 
to fluctuations and the detector resolution not implemented in the model, is not perfectly well 
described. Even replacing the black-disk nucleon-nucleon overlap function with a Gaussian 
does not improve the fit, as the difference in the A^ part distribution is washed out in the 7V C h distri- 
bution. However, it is important to remark that the fit is used solely to determine the AP, which 
is quite insensitive to the detailed shape of the high multiplicity tail. 

An equivalent procedure was applied to fit, with the NBD-Glauber method, the distribution of 
the hits collected in the outer layer of the SPD, and the tracks reconstructed in the TPC. All 
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these analyses give consistent results, which are summarized in Sec. [43 



c\T 9- 




Fig. 11: Centrality dependence of d/V C h/dT] per participant pair as a function of N parU measured in the 
Pb-Pb data at y^NN = 2.76 TeV fitted with various parametrizations of N part and N co \\, calculated with 
the Glauber model. The fit parameters are given in the figure. Data are from ll35l . 

4.2.1 Ancestor dependence 

The number of emitting sources Af ance stors is determined by a function inspired by the two- 
component models, i.e. Af ance stors = / • N par t + (1 — /) • N co \\. However, other assumptions can 
be made leading to a different parametrization, which are briefly discussed in the following. 
The ancestor dependence on /V part and N co \\ derives from a parametrization of the dependence 
of the charged particle multiplicity on /V part and N co \\. Systematic studies of this dependence 
performed at the SPS I55H551. at RHIC 01, and recently at the LHC G30@BD, have been 
used in an attempt to constrain different models of particle production. 

Table 2: Parameters of the fit to the charged particle multiplicity for the three different parametrizations 
discussed in the text, with error and # 2 /NDE 



The charged particle multiplicity is expected to scale with /V part in scenarios dominated by soft 
processes. In this case, all the participant nucleons can be assumed to contribute with the same 
amount of energy to particle production, and the scaling with /V part is approximately linear. 
By contrast, a scaling with N co \\ is expected for nuclear collisions in an energy regime where 
hard processes dominate over soft particle production. In this case, nuclear collisions can be 
considered as a superposition of binary nucleon-nucleon collisions. Two-component models 
are used to quantify the relative importance of soft and hard processes in the particle production 
mechanism at different energies. 



Model 



Normalization Error Fit Par. Error # 2 /NDF 



/•N part + (l-/).N co u 
iv a 

iN| part 
Ncc/ 



2.441 0.281 / = 0.788 0.021 0.347 
1.317 0.116 a = 1.190 0.017 0.182 
4.102 0.297 j8 = 0.803 0.012 0.225 
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To determine the scaling behavior of the particle production, the charged particle multiplicity 
dTVch/dr/ as a function of the number of participants Af par t was fitted with a power-law function 
of Afp ar t i.e. d7V c h/ dr/ oc A^ rt . While at SPS energy the scaling with 7V part is approximately linear, 
i.e. a rsj 1 P6H3811 . results from the experiments at RHIC show evidence of a large contribution 
of hard processes to particle production, resulting in a > 1 . 

The charged particle multiplicity per participant pair dA^h/di] / (0.57V part ) measured by AL- 
ICE [35] is fitted (Fig. [TT]) with three different parametrizations of the ancestor dependence 
mentioned above: 

— a two-components model: d/Vch/drj oc f . Af part + (1 — /) • N co \\; 

— a power-law function of Af part : d/V c h/dr/ oc A^ rt ; 

— a power-law function of N co \\: dA^/dr/ oc A/^ oll . 

and the fit parameters are reported in Tab. [2} We note that the value obtained for / is in a good 



agreement with the value obtained in the NBD-Glauber fit, shown in Fig. 10 



While the value obtained for a and for j8 with the power-law parametrization of 7V part and 
A/coii indicate that neither of these scalings perfectly describes the data (a > 1 and j8 < 1), we 
note that the value of a is similar to that measured at RHIC (1.16 ± 0.04 1(391 ) and slightly 
higher than that at the SPS (a ~ 1, see ll36l for a review). The results obtained with the two- 
component model, where < / < 1, indicate that both the contribution of 7V part and N co \\ are 
needed to explain the particle production confirm this. However, the # 2 /NDF reported in Table 
[2j indicate an equally good fit for all models, thus revealing that no unique physics conclusion 
can be drawn from such fits and that the particular choice of parametrization has no influence 
on the results of the centrality determination. 

4.3 Systematic uncertainty on the Anchor Point 

The determination of the AP by either correcting or fitting the multiplicity distribution is eval- 
uated in Fig.[12]by comparing the VZERO amplitude distributions for various event selections. 
The systematic uncertainty on the AP is estimated by comparing the percentage of the hadronic 
cross section at the VZERO amplitude chosen as the AP (VOap) obtained correcting or fitting the 
multiplicity distribution. For the first method (Sec. |4.1[ ), we used the results from the HIJING 



and AMPT simulations. For the second method (Sec. |4.2| ), we used alternative centrality defi- 
nitions based on (i) TPC tracks; (ii) SPD hit multiplicities, and obtaining a value for the VOap 
using the correlation between SPD or TPC and VZERO; (iii) different ranges for the Glauber 
model fit; (iv) different ancestor dependence of the particle production model to a power law of 
Af par t; (v) different nucleon-nucleon cross section and parameters of the Woods-Saxon distribu- 
tion within their estimated uncertainties. All the results, compared in Table[3j allow to define the 
AP as the VZERO amplitude above which we obtain 90% of the hadronic cross section with the 
NBD-Glauber fit (the baseline in Tabj3]) with a systematic uncertainty of 1%, determined as the 
RMS of all the results presented in Table [3j The variations of the AP are not part of the quoted 
systematic uncertainties for Af part , N co \\ and Taa, which include only variations of the Glauber 
parameters. The uncertainty on the AP is typically included in our analyses as an uncertainty on 
the limits of the centrality classes and propagated into an uncertainty on the specific measured 
observable. 
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Fig. 12: Distribution of the VZERO amplitude zoomed in the most peripheral region. The distribution is 
compared to the NBD-Glauber fit and to the sum of the HIJING + STARLIGHT + QED simulations. 

5 Centrality classes and their relation with geometrical quantities 

5.1 Determination of the centrality classes with the multiplicity distributions 

The percentile of the hadronic cross section is determined for any value of the VZERO am- 
plitude by integrating the measured VZERO amplitude distribution normalized at the anchor 
point VOap, i.e. 90% of the hadronic cross section. For example, if we define V as the VZERO 
amplitude, the top 10% central class is defined by the boundary VOio which satisfies 

Jvo w (dN en /dV)dV = l 
IvoJdN evt /dV)dV 9 

The same is done for the number of clusters in the SPD and the number of reconstructed tracks 
in the TPC. The events with multiplicity lower than that of the anchor point, contaminated by 
EM background and trigger inefficiency, are not used in the physics analyses. 

One can divide the experimental distribution into classes, by defining sharp cuts on e.g. VZERO 
amplitude, which correspond to well defined percentile intervals of the hadronic cross section. 
The number of centrality classes that one can define is connected with the resolution achieved 
on the quantities used in the definition. In general, centrality classes are defined so that the 
separation between the central values of b and 7V pa rt for two adjacent classes is significantly 
larger than the resolution of that variable (see Sec. [6]). 

5.2 Finding the number of participants with the multiplicity distributions 

In Sec. |4.2| we fit the measured VZERO amplitude distribution with the amplitude distribution 
simulated with the NBD-Glauber. This creates a connection between an experimental observ- 
able and the geometrical model of nuclear collisions used in the Glauber Monte Carlo. From 
this we can access the geometrical properties, like A^ ar t, N co \\, Taa- A given centrality class, de- 
fined by sharp cuts in the measured distribution, corresponds to the same class in the simulated 
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Table 3: Comparison of the percentage of the hadronic cross section above the VZERO amplitude chosen 
as AP (VOap) for various cases considered in the systematic studies of the Glauber fits and with HIJING 
and AMPT simulations. 

Method % of total cross section above the VOap 

Glauber Fits 
Baseline 90.00 



(i) fit TPC tracks 


89.88 


(ii) fit SPD clusters 


89.87 


(iii) fit only 50% of cross sec 


90.11 


(iv) different ancestor dependence 


90.66 


(v) different Wood-Saxons par 


90.43 


HIJING simulations 




2— out — of— 3 


92.50 


VOAND 


89.05 


3— out — of— 3 


90.15 


VOAND + TPC 


91.12 


VOAND + ZDC 


89.52 


AMPT simulations 




2— out — of— 3 


92.49 


VOAND 


89.49 


3— out — of— 3 


90.59 


VOAND + TPC 


91.36 


VOAND + ZDC 


89.00 



distribution. For the simulated distribution we retain the input information from the Glauber 
model. Therefore, we can calculate the mean number of participants (Af par t) 9 the mean number 
of collisions (N co \\) and the average nuclear overlap function (7aa) f° r centrality classes de- 
fined by sharp cuts in the simulated multiplicity distribution, corresponding to given percentiles 



of the hadronic cross section. As shown in Table |A.1[ the mean values and their dispersions 
differ from those calculated for geometrical classes, defined by sharp cuts in the impact param- 
eter b (Table [TJ, by less than 1% for the most central classes (up to about 50%) and by less than 
2% for the most peripheral ones (above 50%). This confirms that multiplicity fluctuations and 
detector resolution only play a minor role in the centrality determination. 

5.3 Determination of the centrality classes with the ZDC 

Another way to determine the centrality is to measure the energy deposited by the spectators in 
the ZDC. The spectator neutrons and protons having a rapidity close to that of the beam, are 
detected in the ZDC. Naively, measurement of the number of spectator neutrons and protons 
would give direct measurement of the number of participants since Af part would simply be given 
by 

N paYt = 2A-E ZDC /E A (7) 

where Ezdc is the energy measured in the ZDC, A = 208 is the mass number of Pb, and 
is the beam energy per nucleon. However, fragment formation amongst the spectator nucleons 
breaks the simple linear and monotonic relation in the measured variables, since some spectator 
nucleons are bound into light nuclear fragments that have a charge over mass ratio similar to 
the beam, therefore, remaining inside the beam-pipe and are undetected by the ZDC E2ll43ll . 
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ZEM amplitude (a.u.) 



Fig. 13: Spectator energy deposited in the ZDC calorimeters as a function of ZEM amplitude. The 
same correlation is shown for different centrality classes (5%, 10%, 20% and 30%) obtained by selecting 
specific VZERO amplitudes. The lines are a fit to the boundaries of the centrality classes with linear 
functions, where only the slope is fitted and the offset point is fixed (see text). 

This effect becomes quantitatively important for peripheral events and therefore Eq. [7] cannot 
be used as a reliable estimate of Af part . 

Consequently, the ZDC information needs to be correlated to another quantity that has a mono- 
tonic relation with 7V part . In our case, we use the energy measured by two small EM calorimeters 
(ZEM). These detectors are placed only on the A side about 7.5 m from the interaction point, 
covering the region 4.8 < r/ < 5.7 ifTOll . 

Since the ZDC calorimeters are far from the interaction region, and therefore have an acceptance 
insensitive to the vertex position, a centrality measurement based on the ZDC is particularly 
suited for any analysis that does not require a vertex cut ll44ll . 

Centrality classes are defined by cuts on the two-dimensional distribution of the ZDC energy 
as a function of the ZEM amplitude. The ZDC signal is proportional to 7V part for central events, 
while the ZEM amplitude is an unknown function of 7V part and N co \\. Therefore the definition 



of the centrality classes in this two-dimensional space is not trivial. As shown in Fig. 13 
centrality classes are defined by using the centrality classes defined previously with the VZERO 
amplitude to determine regions in the ZDC-ZEM plane, corresponding to a given centrality. The 
boundaries between centrality classes, or the points belonging to the same narrow centrality 
class c (c ± 8c) can be fitted with linear functions. All these lines are found to intersect at a 
common point. Using this common point, we refitted the boundaries of the various centrality 



classes with the linear functions shown in Fig. 13 



As can be seen from the figure, the slopes of the fitted functions increase going from central 
to peripheral collisions and tend to infinity, as the lines become almost straight vertical lines, 
when approaching the point where the correlation between ZDC and ZEM inverts its sign. The 
value of the slope that defines a centrality class in the ZDC vs ZEM phase-space is proportional 
to the tangent of the percentile, which implies that the percentiles behave like an angle in the 
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5000 10000 15000 20000 25000 

VZERO a mplitude (a.u.) 



Fig. 14: VZERO amplitude distribution of events of various centrality classes selected from the correla- 
tion between ZDC and ZEM amplitudes explained in the text. 

ZDC vs ZEM phase-space. 

This function of ZDC and ZEM can then be used as centrality estimator for the most central 
events (0-30%) above the turning-point of ZDC. Figure[14]shows the distribution of the VZERO 
amplitude for all triggered events and for various centrality classes selected with this method. 



6 Resolution of the centrality determination 

As described above, two independent methods are used to determine experimentally the cen- 
trality of the collision. The first one uses the multiplicity distributions from various detectors 
covering different pseudo-rapidity ranges. Specifically we use the sum of the amplitude in the 
VZERO detectors (A and C side), the number of clusters in the outer layer of the SPD detector, 
and the number of tracks reconstructed in the TPC. The second method uses the ZDC correlated 
with the ZEM. 

The accuracy of the experimental determination of the centrality was evaluated by comparing 
the different estimates event-by-event. For example, in Fig. [15] we compare the estimates based 
on the SPD multiplicity and the VZERO amplitude. The VZERO amplitude distribution is 
shown for two centrality classes selected by the SPD multiplicity. The distributions for the two 
centrality classes are reasonably well fitted with a Gaussian distribution. 

The resolution in the experimental definition of the centrality classes is evaluated event-by- 
event as the RMS of the distribution of the differences between the centrality determined over 
all estimators and the mean value of the centrality for the event. First the average value of 
the centrality (c) is calculated for each event by averaging the centrality determined by each 
estimator: 

(c) = Sk* (8) 

C[ is the centrality of an event determined by an estimator i 9 where i is the index running over 
all N = 6 centrality estimators used: VZERO (A and C), SPD, TPC, ZDC. In the next step, the 
centrality is weighted by A/ = q — (c) : the difference between the centrality determined by each 
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VZERO amplitude (a.u.) 



Fig. 15: Top: Correlation between SPD multiplicity and VZERO amplitude. The rapidity coverage of 
each detector is indicated on the figure. Bottom: VZERO amplitude distributions for the centrality classes 
selected by SPD. Two centrality classes (1-2% and 40-45%) are indicated and fitted with a Gaussian. 



estimator and the mean value of the centrality from Eq. [8j 



(9) 



This latter calculation is performed iteratively replacing (c) by the new value until convergence 
is achieved which typically occurs after the second iteration. Finally, the centrality resolution 
of an estimator is evaluated as the RMS of its A/ distribution for each centrality. 

The ZDC-ZEM estimator is ignored for peripheral events ((c) > 35%) since its results are re- 



liable only for the most central collisions. The resolution is shown in Fig. [16] (left panel) as a 
function of the centrality percentile. 

The resolution depends on the rapidity coverage of the detector used. The best centrality res- 
olution is achieved when combining the VZERO-A and VZERO-C detector, due to the large 
pseudo-rapidity coverage (4.3 units in total). It ranges from 0.5% in central to 2% in peripheral 
collisions. The resolution obtained with the SPD and the TPC ranges from 1% in central to 3% 
in peripheral collisions (> 80%). 

We measured the pseudo-rapidity dependence of the charged particle multiplicity at midrapid- 
ity Il45ll with the SPD, and at forward rapidity [46] using all the rapidity coverage of the SPD, 
the VZERO and the FMD detectors. The total charged particle multiplicity is obtained by 
integration. The centrality resolution was scaled by y/N C h measured in the rapidity window of 



each detector (see right panel of Fig. [16]). The figure shows that all the results are consistent on 
an arbitrary unit scale, except for the ZDC-ZEM estimator which is better for central collisions 
because it uses information from two detectors. 
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Fig. 16: Left: Centrality resolution A z for all the estimators evaluated in the analysis. Right: Resolution, 
in arbitrary units, scaled by y/N^ measured in each detector. 

The centrality resolution was tested with a full HIJING and GEANT detector simulation. In 
the HIJING simulations the true value of the event centrality (c trU e) is known for every given 
event. After using GEANT one obtains the signals in VZERO, SPD, TPC for the given event 
and hence using these centrality estimators can calculate the value of the (c) for the given event 
with Eq|8j The real centrality resolution, given for the given event by the difference between the 
Ctrue and the (c) calculated for each extimator, is consistent with the one calculated with data. 



7 Summary 

Heavy-ion collisions can be characterized by the number of charged particles produced in the 
collision. In principle, when normalized to the trigger efficiency used to collect the data sample, 
the charged particle multiplicity could provide a measurement of the hadronic cross section. 
However, at the LHC the large cross section for EM processes contaminates the very peripheral 
collisions. This problem was overcome in two ways. 

In the first method, dedicated simulations of hadronic and EM processes (Fig. [8]) were per- 
formed and data were corrected for efficiency of the event selection (Fig. [7]) and purity of the 
event sample (Fig. [9]). In the second method, the measured multiplicity distribution was fitted 



with a Glauber calculation (Fig. [TO]). Both methods allow to determine a centrality value above 
which the background contamination is negligible and the event selection is fully efficient. The 
corresponding value of multiplicity and centrality is defined as the anchor point, and is used for 
the centrality normalization (Table [3]). 

Using the AP, the measured event sample can be divided in centrality classes which correspond 
to well defined percentiles of the hadronic cross section. Several approaches were developed. 
The first method uses charged particle multiplicity (measured by various detectors, with differ- 
ent rapidity coverage, such as the VZERO, the SPD, and the TPC). The second method uses 
the ZDC, which measures the nucleon spectators directly, as well as the correlation to the ZEM 
energy in order to resolve the ambiguity due to nuclear fragmentation. The centrality is ob- 
tained from linear functions that fit the contours of the classes defined by the VZERO, in the 



ZDC-ZEM plane (Fig. [13]). As standard method, typically used in ALICE physics analyses, we 
used the NBD-Glauber fit to the VZERO amplitude (Fig. [TO]) to determine the AP, and the other 
methods described to asses a systematic uncertainty on the centrality determination. 
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The resolution of the centrality determination, which depends on the pseudo-rapidity coverage 
of the detector used, was determined as the weighed RMS of all the estimates; it ranges from 



0.5% in central to 2% in peripheral collisions (Fig. 16). 



Finally, mean numbers of the relevant geometrical quantities, such as Af part and N co \\, were 
calculated for typical centrality classes, using the Glauber Model and the fit to the measured 
multiplicity distribution (Table [T]). This fit creates a mapping between a measured quantity 
and one obtained with a phenomenological calculation for which the geometrical properties 
are known. The results, nearly identical to those obtained for centrality classes defined by 
classifying the events according to their impact parameter, provide a general tool to compare 
ALICE measurements with those of other experiments, at different energies and with different 
colliding systems as well as theoretical calculations. 
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A Tables 

As described in Sec. [2j for the physics in ALICE analyses the average values of Af part , N co \\, 
or Taa for centrality classes defined by sharp cuts in the impact parameter distributions are 
used. These are reported in Table [TJ Therefore (Af part ), (N co \\) an d (Taa) depend exclusively 
on the nuclear geometrical parameters, and not on any measured quantity. Their uncertainty 
is calculated by varying the parameters of the Glauber calculations (i.e. the parameters of the 
Woods-Saxon and the hadronic cross section Oj^ 1 ) by the known uncertainty. We label the 
(TVpart) calculated with this procedure as (A^°). 

Another possibility, discussed in Sec ]5.2[ is to define the average values of Af part , A/coii, or Taa 
for centrality classes by sharp cuts in the fitted multiplicity distribution. Following this strat- 
egy, it is also possible to incorporate in the uncertainty, besides the uncertainties related to the 
Glauber calculation, those related to the measurement of the AP: the experimental region which 
is actually being used for the physics analyses, because it is free of background and the trigger 
efficiency is known. In this case, the AP can be varied by the uncertainty that was estimated 
(90% ±1%) and recalculate (Af part ), (N co \i) and (Taa) with these variations. The (Af part ) calcu- 
lated with this procedure is labled as (N^). The variations for the AP are labelled (A^^ t a+ ) 
and (Afpart - ) respectively. 

(A^ a e r ° t ) is compared to (N^) for various centrality classes. The default values 



A.l 



In Table 

of (NpHfj are compared to the values obtained by varying the AP. The discrepancies A are 
calculated as 



Same comparison is done in Tables A.2 and A.3|for (7V co n)and (rAA)respectively. Table A.4 



gives the comparison for the three quantities but for bigger centrality classes. 
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Table A.l: Af par t for Pb-Pb collisions at y^NN = 2.76 TeV with the corresponding uncertainties derived 
from a Glauber calculation. The (A^f) are calculated from the NBD-Glauber fit to the VZERO ampli- 
tude, while the (A^ r °) are obtained by slicing the impact parameter distribution. (A^f) is also calculated 
for two variations of the AP, i.e. moving it to 91% ((A^ r \ a+ )) and to 89% ((A^ar t a+ )) respectively. The last 
three columns report the discrepancies between (A^ r °) and (A^f) and (Afp^ t a )with the uncertainty of the 
AP. 
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35-40% 


116(2.2) 


116.7 


11 


118.6 


114.7 


0.31 


0.81 


0.86 


40-45% 


94.11 (2.1) 


94.77 


9.7 


96.68 


92.83 


0.35 


1 


1 


45-50% 


75.3 (2.3) 


75.91 


8.4 


77.72 


74.02 


0.4 


1.2 


1.3 


50-55% 


59.24 (2.5) 


59.77 


7.3 


61.49 


58.02 


0.44 


1.4 


1.5 


55-60% 


45.58 (2.9) 


46.1 


6.3 


47.66 


44.47 


0.57 


1.7 


1.8 


60-65% 


34.33 (2.6) 


34.65 


5.4 


36.09 


33.2 


0.47 


2 


2.1 


65-70% 


25.21 (4) 


25.38 


4.5 


26.62 


24.16 


0.34 


2.4 


2.5 


70-75% 


17.96 (3.3) 


18.06 


3.8 


19.07 


17 


0.27 


2.7 


3 


75-80% 


12.58 (3.7) 


12.45 


3 


13.25 


11.61 


0.54 


3.1 


3.5 


80-85% 


8.812(2.8) 


8.275 


2.4 


8.914 


7.646 


3.1 


3.7 


4 


85-90% 


6.158(2.4) 


5.516 


1.8 


6.035 


3.406 


5.5 


4.5 
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Table A.2: Same as Table |Al] for A^ u. 



Cent. 


<AQ (syst. %) 




RMS 






(%) 


(%) 


(%) 


0-1% 


1861 (8.2) 


1863 


83 


1864 


1863 


0.059 


0.016 


0.018 


1-2% 


1766 (8.2) 


1761 


79 


1762 


1759 


0.15 


0.04 


0.047 


2-3% 


1678 (8.2) 


1678 


79 


1680 


1676 


0.0063 


0.063 


0.067 


3-4% 


1597 (8.3) 


1596 


78 


1599 


1592 


0.039 


0.096 


0.1 


4-5% 


1520 (8.1) 


1520 


77 


1524 


1516 


0.002 


0.12 


0.12 


5-10% 


1316(8.2) 


1316 


110 


1321 


1310 


0.0099 


0.2 


0.2 


10-15% 


1032 (8.2) 


1034 


95 


1040 


1027 


0.083 


0.32 


0.34 


15-20% 


809.8 (8) 


811.7 


80 


819.2 


803.9 


0.11 


0.46 


0.48 


20-25% 


629.6 (7.8) 


631 


68 


639 


622.9 


0.11 


0.63 


0.65 


25-30% 


483.7 (7.5) 


485.8 


57 


493.6 


477.7 


0.22 


0.79 


0.85 


30-35% 


366.7 (7.4) 


368.4 


49 


375.9 


360.7 


0.23 


1 


1.1 


35-40% 


273.4 (7.4) 


274.8 


40 


281.7 


267.8 


0.26 


1.2 


1.3 


40-45% 


199.4 (6.9) 


200.7 


33 


206.8 


194.5 


0.32 


1.5 


1.6 


45-50% 


143.1 (6.6) 


143.8 


26 


149 


138.6 


0.26 


1.8 


1.9 


50-55% 


100.1 (6.5) 


100.6 


20 


104.9 


96.22 


0.25 


2.1 


2.2 


55-60% 


68.46 (6.2) 


68.7 


15 


72.12 


65.28 


0.18 


2.4 


2.6 


60-65% 


45.79 (5.7) 


45.79 


12 


48.47 


43.12 


0.0038 


2.8 


3 


65-70% 


29.92 (6.8) 


29.66 


8.3 


31.68 


27.73 


0.43 


3.3 


3.4 


70-75% 


19.08 (5.7) 


18.82 


5.8 


20.22 


17.39 


0.68 


3.6 


4 


75-80% 


12.07 (5.7) 


11.62 


4 


12.61 


10.63 


1.9 


4.1 


4.4 


80-85% 


7.682 (5.1) 


6.925 


2.7 


7.595 


6.269 


5.2 


4.6 


5 


85-90% 


4.904 (4.1) 


4.148 


1.8 


4.651 


2.257 


8.4 


5.7 
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Table A.3: Same as Table [AT] for T AA . 



Cent. 


/T geo \ (svst %) 
x^AB / W» L - 10 ) 


/T^datax 
X^AB / 


RMS 


/7^data+\ 
^AB / 


/ 7^ data— \ 
X^AB / 


a data (C7\ 


Kdata+ (c?\ 
^data V /c ) 


sdata- rci\ 
^data V /c ) 


u-i /c 


9Q 08 (^K 1\ 


9Q 1 1 
Ly . 1 1 


i .j 


9Q 1 9 
Ay. 1Z 


9Q 1 
Ly . 1 


0^£ 
U.UJO 


U.UIO 


018 
U.UIO 


1 9% 
i-z /o 


97 £ 9"> 
Z / .o ^J.zj 


97 ^1 
Z/.Jl 


1 9 
1 .Z 


97 ^ 

Z / Jj 


97 ZL8 
Z / .^fo 


O 1 f\ 


OZL 


0ZL7 


9 

ZO /C 


9^ 99 c\ 9\ 

ZD.ZZ ^J.Z ^ 


9£ 99 

ZD.ZZ 


1 9 
1 .z 


9£ 9^ 
ZD.ZJ 


9£ 1 Q 
ZD. 


00^8 


C\f\^ 
U.UDJ 


0£»7 
U.UD / 


J -4- /C 


9ZL CK 9"i 


94 en 


1 9 
1 .z 


9ZL Q8 


9ZL 88 


0^9 


OQ^» 


1 

U. 1 


H— J /C 


9a 7< /a 9\ 
ZJ. / J V-j.z^ 


9^ 7S 
ZJ. / J 


1 9 
1 .z 


9^ 8 


9^ £»Q 


001 Q 


1 9 
U.1Z 


1 9 
U. 1Z 


s 10% 


90 C\ 'W 

ZU.JO ^J..J ^ 


90 


1 8 
1 .o 


90 f\A 


90 ZL7 


00^6 


9 
u.z 


9 
u.z 


iu-ij /c 


ID. 1 J \J Sj) 


1 £ 1 ^ 
ID. 1 J 


i < 
1 .J 


1 fx Ifx 

lO.ZO 


1 £ OzL 


0^»8 
U.UOo 


O ^9 
U. JZ 




1 ^ 90% 


19 n 7\ 


1 9 £8 


1 9 
1 .Z 


1 9 8 
IZ.o 


1 9 

IZ.JD 


O 1 ^ 


Af\ 


O ZL8 


90 9^% 


7.0J / ^J. / ) 


Q 8£» 
7.0D 


1 1 
1 . 1 


Q Q84. 


Z?. / J J 


1 9 
U. 1Z 


£\"k 
U.DJ 


£»S 
U.DJ 


9< ^0% 

ZJ- D\J /C 


7 SS8 A\ 


7 ^Q1 

/ .Jy 1 


Q 


7 71^ 
/ . / 1 J 


/ .H-D J 


99 
u.zz 


7Q 
U. / y 


8S 


^0 


S 7^ n ^ 

J. / D \D. J) 


S 7S6 

J . / JO 


lf\ 


S 87^ 


J .OJO 


9^ 

U.Zj 


1 


1 1 

1 . 1 


35-40% 


4212 (3.7) 


4.294 


0.63 


4.402 


4.184 


0.26 


1.2 


1.3 


40-45% 


3.115(3.9) 


3.136 


0.51 


3.231 


3.039 


0.33 


1.5 


1.6 


45-50% 


2.235 (4.2) 


2.248 


0.41 


2.328 


2.165 


0.28 


1.8 


1.9 


50-55% 


1.564 (4.7) 


1.572 


0.32 


1.639 


1.504 


0.25 


2.1 


2.2 


55-60% 


1.07 (5.2) 


1.073 


0.24 


1.127 


1.02 


0.16 


2.4 


2.6 


60-65% 


0.7154 (5) 


0.7154 


0.18 


0.7573 


0.6737 


0.0007 


2.8 


3 


65-70% 


0.4674 (6.2) 


0.4635 


0.13 


0.4949 


0.4333 


0.42 


3.3 


3.4 


70-75% 


0.2981 (6.4) 


0.2941 


0.091 


0.3159 


0.2717 


0.68 


3.6 


4 


75-80% 


0.1885 (6.9) 


0.1815 


0.062 


0.197 


0.1661 


1.9 


4.1 


4.4 


80-85% 


0.12(6.5) 


0.1082 


0.042 


0.1187 


0.09795 


5.2 


4.6 


5 


85-90% 


0.07662 (5.9) 


0.06481 


0.028 


0.07267 


0.03526 


8.4 


5.7 
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Table A.4: Same as above with bigger centrality classes 



Cent. 


(N^) (syst. %) 




RMS 


«t a+ > 


«r> 


A data (%) 






0-5% 


382.7 (0.77) 


382 


17 


382.3 


381.7 


0.096 


0.04 


0.038 


5-10% 


329.6 (1.3) 


329.7 


18 


330.5 


328.8 


0.015 


0.13 


0.13 


10-20% 


260.1 (1.5) 


260.5 


27 


261.9 


259.1 


0.079 


0.26 


0.27 


20-40% 


157.2 (2) 


157.8 


35 


159.6 


155.9 


0.19 


0.58 


0.6 


40-60% 


68.56 (2.9) 


69.13 


22 


70.89 


67.35 


0.42 


1.3 


1.3 


60-80% 


22.52 (3.4) 


22.64 


12 


23.76 


21.51 


0.27 


2.4 


2.6 




(A© (syst. %) 




RMS 


« + ) 


«-) 


A data (%) 


Agf (%) 


K7a~ (*) 


0-5% 


1685 (11) 


1684 


1.4e+02 


1686 


1681 


0.044 


0.065 


0.067 


5-10% 


1316(11) 


1316 


l.le+02 


1321 


1310 


0.011 


0.2 


0.2 


10-20% 


921.2(10) 


922.7 


1.4e+02 


929.8 


915.3 


0.079 


0.39 


0.4 


20-40% 


438.4 (9.7) 


440 


1.5e+02 


447.5 


432.3 


0.19 


0.85 


0.89 


40-60% 


127.7 (8.8) 


128.4 


59 


133.2 


123.7 


0.29 


1.8 


1.9 


60-80% 


26.71 (7.3) 


26.48 


18 


28.25 


24.74 


0.43 


3.2 


3.4 




(7^°) (syst. %) 


/xdata\ 
^AB / 


RMS 


/7^data+\ 
X^AB / 


/T^data— \ 
X^AB / 


A data (%) 




K7a~ (*) 


0-5% 


26.32 (3.2) 


26.31 


2.2 


26.34 


26.27 


0.028 


0.066 


0.066 


5-10% 


20.56 (3.3) 


20.56 


1.7 


20.64 


20.47 


0.0051 


0.2 


0.2 


10-20% 


14.39 (3.1) 


14.42 


2.2 


14.53 


14.3 


0.092 


0.39 


0.4 


20-40% 


6.85 (3.3) 


6.876 


2.3 


6.993 


6.754 


0.19 


0.85 


0.89 


40-60% 


1.996 (4.9) 


2.007 


0.92 


2.081 


1.933 


0.28 


1.8 


1.9 


60-80% 


0.4174 (6.3) 


0.4137 


0.29 


0.4414 


0.3865 


0.44 


3.2 


3.4 
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